Evidence for ±s-wave pairing symmetry in LiFeAs: specific heat study 
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We report specific heat capacity measurements on a LiFeAs single crystal at temperatures down 
to 400 mK and magnetic fields up to 9 Tesla. A small specific heat jump at Tc and finite residual 
density of states at T = K in the superconducting (SC) state indicate that there are strong 
unitary scatterers that lead to states within the SC gap. A sub-linear magnetic field dependence of 
the Sommerfeld coefficient 'y{H) at T = K is equally well fitted by both a nodal d-wave gap as well 
as a sign changing multiband ±s-wave gap. When impurity effects are taken into account, however, 
the linear temperature dependence of the electronic specific heat Cd/T at low temperatures argues 
against a nodal d-wave superconducting gap. We conclude that the SC state of LiFeAs is most 
compatible with the multiband ±s-wave SC state with the gap values As^aj; = 0.46Aiarae- 

PACS numbers: 74.70.Xa, 74.25.Bt, 74.20.Rp 



I. INTRODUCTION 



Information on the superconducting (SC) pairing sym- 
metry is important for the understanding of the SC pair- 
ing mechanism. For the newly discovered Fe-pnictide 
superconductors, determination of the SC order param- 
eter has been controversial partly due to its sensitive 
dependence on measurement probes and material stoi- 
chiometryivS, in the '1111' phase of REFeAsO (RE = 
rare earth), a fulUy gapped scenario is supported by 
angle-resolved photoemission spectroscopy(ARPES) and 
tunnelling measurements)^!^ but a gap symmetry with 
nodes on the Fermi surface is suggested by far-infrared el- 
lipsometric measurements.— Penetration depth measure- 
ments reported either a full gap or a nodal gap for RE 
= Sm,Pr&- and RE=Laj^ respectively. Identification of 
the SC gap symmetry is also controversial in the '122' 
phase of AeFe2As2 (Ae = alkaline earth). For electron 
doped BaFe2As2, full gaps are reported from ARPES^ 
and scanning tunnelling microscopy (STM) measure- 
mentSj^- but anisotropic or nodal gaps are supported 
by specific heat,— thermal conductivity ) ^^i "'^^ penetration 
depth, ^'14 Raman scatteringfi^ nuclear magnetic reso- 
nance (NMR),-i^ and muon spin resonance (/J,SR).^- For 
hole-doped Bal22, multiple gaps without nodes are re- 
ported from specific heaii^ and thermal conductivity.— 

LiFeAs, the so-called '111' phase, is unique in that it is 
superconducting without carrier doping and the residual 
resistivity ratio (RRR) is about 50, one of the highest 
among Fe-based SC compounds,^" holding promise for 
a determination of the superconducting gap symmetry 
without complications from material defectSi ^^'^^ Various 
spectroscopic and thermodynamic measurements suggest 
an isotropic SC gap without nodes on the Fermi sur- 
facei^"— Previous specific heat measurements, a direct 



bulk probe of the electronic density of states (DOS), sug- 
gest multiple SC gaps to explain its temperature depen- 
dence!^^ The analysis, however, is limited due to lack 
of specific heat data below 2 K, where a difference in 
the DOS becomes prominant among different types of 
SC pairing symmetry. In order to elucidate the nature of 
the order parameter of LiFeAs, we measured the specific 
heat of LiFeAs single crystals down to 400 mK and under 
magnetic field up to 9 Tesla. A large value of the Som- 
merfeld coefficient 70 (~ 7.5 mJ/mol-K^) in the SC state 
at zero magnetic field and small specific heat jump ratio 
AC/C„ at Tc{~ 0.5) indicates the presence of large elas- 
tic scattering due to impurities. Here C„ is the specific 
heat at Tc in the normal state. When the impurity effects 
are taken into account, observation of a linear tempera- 
ture dependence of the low-T specific heat C /T and a 
sublinear magnetic field dependence of the Sommerfeld 
coefficient ^{H) excludes a d-wave SC gap, but is con- 
sistent with the sign changing multiband ±s-wave state 

with the gap values Asmall = OAGAlarge- 



II. EXPERIMENTS 

Single crystalline LiFeAs, which is formed in a 
P4/nmm tetragonal Cu2Sb-type structure with a — 
3.7818 A and c = 6.3463 A, was synthesized in a sealed 
tungsten crucible by the Bridgeman method;^ X-ray 
diffraction pattern analysis showed that crystals from 
this batch are homogeneous and well oriented. This 
growth technique seems to avoid inclusion of impurities 
that can lead to an anomalous Schottky-like upturn in 
the low-temperature specific heat of Fe-based supercon- 
ductors, making it a useful synthesis technique for a spe- 
cific heat study of SC gap symmetry.— The crystals used 
in the current study are from the same batch that has 
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FIG. 1: (color online) (a) Specific heat of LiFeAs at zero mag- 
netic field is plotted down to 400 mK. The red solid line is an 
estimated normal state specific heat Cn/T from a polynomial 
form described in the text, (b) Calculated electronic entropy 
in the superconducting state (Ses) and in the normal state 



been used for previous transport, thermodynamic, and 
spectroscopic studies i^°i^^'^^ A Quantum Design PPMS 
(Physical Properties Measurement System) with a '^He 
option was used to measure specific heat down to 400 mK 
and up to 9 T for magnetic fields applied parallel and 
perpendicular to the crystalline c-axis. 



III. RESULTS 

Figure 1(a) shows the temperature-dependent specific 
heat capacity of a LiFeAs single crystal at zero field. A 
superconducting transition occurs at 15 K {=Tc), which 
is determined from the mid-point of the the specific heat 
anomaly. The measured specific heat (C) consists of 
electronic {Cei) and phononic (Cph) contributions, i.e., 
C = Cei + Cph- Non-Debye behavior has been often re- 
ported in Fe-based superconductors and ascribed to a 
large Einstein contribution or a low Debye temperature 
Td. Because an Einstein phonon contribution is negligi- 
ble below 25 K in LiFeAs, a second term of the harmonic- 
lattice approximation is added to the phonon specific 
heat: C„ = 7n(0)T + PiT^ + /32r'5,27 The specific heat 
is best described by the least-squares fit with 7„(0) — 
17.0 ± 0.9 mJ/mole-K^, /3i = 0.231 ± 0.004 mJ/mole-K^ 
and P2 — —0.0001 mJ/mole-K^. The electronic Sommer- 
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FIG. 2: (color online) Normalized electronic specific heat co- 
efficient (black open circles) and the best fits of the two SC 
gap models: d-wave (blue circles) and ±s-wave (red trian- 
gles) models. For the d-wave model, 2l\rnax/Tc = 3.5 and 



r = Wi, 



o/'^^^norj 



— 0.13 are obtained for the best re- 



sults. Here, Amax is the maximum gap of the d-wave gap 
A(6I). For the ±s-wave model, |A|smaii/|A|iar9e = 0.46, 



2A,a,.5,/rc = 3.8, and F 



o/ivNtotai = 0.17 are obtained 



from the least squares fit. We note that the inversely related 
ratios between \A\s,naii/\A\iargB and N small /Niarge are an in- 
trinsic property of the ±s-wave pairing model with a domi- 
nant interband paring interaction!^ Inset: Low- temperature 
specific heat in the main panel is magnified to show a notable 
difference between the two gap models. 



feld coefficient 7„(0)(= 17.0 ± 0.9 mJ/mole-K^) at zero 
magnetic field in the normal state is comparable with 
previously reported values of 23, 20, and 10 mJ/mole-K^ 
from Refs. [28], [24], and [27], respectively. The De- 
bye temperature estimated from the fit is 294 K, which 
is similar to 310 K from Ref. [27]. Figure lb shows 
the electronic entropy in the normal (squares) and su- 
perconducting (circles) states after subtracting the non- 
electronic contribution. The two entropies become equal 
at Tc, which satisfies the entropy constraint, = L" {C— 
Cn)/TdT. From the electronic entropy in Fig. 1(b), 
we estimate the superconducting condensation energy 
U = J^^iSen - Ses)dT = Bc{Qf/2^io -363 mJ/mole and 
thermodynamic critical field -Bc(O) ~0.23 T. When we 
use the Ginzburg-Landau parameter n = Xab/S,ab =29±7 
of LiFeAs from a SANS experiment, ^"^ the upper critical 
field Moi^c2 (= V2kBc{0)) is 9.4±2.3 T. For comparison, 
^qH^2 from penetration measurements is 17 Ti^ 

Figure 2 shows the electronic specific heat coefficient 
of LiFeAs at zero magnetic field after subtracting the 
phononic contribution: jei{T) — Cei/T—{C — Cph)/T. 
The specific heat jump AC at Tc is 7.65 mJ/mol-K^ at 
zero field. The jump relative to the normal state specific 
heat Cn is 0.5 {—AC/Cn), which is much smaller than 
the BCS value 1.43, indicating that the quasiparticles 
participating in the SC condensation experience strong 
elastic scattering because inelastic scattering usually en- 
hances the jump ratio.^'^ Also the substantial value of the 



Sommerfeld coefficient 70 in the SC state, which accounts 
for about 45 % of the normal state value 7„(0), also sup- 
ports the presence of strong (unitary) scatterers because 
weak (Born limit) scatterers are not capable of inducing 
states inside the SC gap of unconventional superconduc- 
tors;^ These results suggest that any model to explain 
the nature of the SC gap symmetry of LiFeAs should take 
into account the effects of impurity scattering. 

To extract information on the gap symmetry from the 
specific heat of LiFeAs, we consider two typical uncon- 
ventional SC gap models that represent two extreme lim- 
its, i.e., a d-wave SC gap with nodes and a ±s-wave gap 
without nodes, and fit the experimental data in Fig. 2. 
Here we used the d-wave model to represent the generic 
behavior of a nodal SC gap state. The impurity effect 
is calculated using the self-consistent 7~-matrix approxi- 
mation (SCTA)."^*^ In each gap model, we determine the 
fitting parameters such as 2A/T'c, the impurity concen- 



tration parameter F 



,/ttN{Q) by finding the best 



overall fit to the data, where A^(0) is the zero-energy den- 
sity of states. The results are shown in Fig. 2 overlayed 
with the experimental data. Both models provided rea- 
sonably good fittings for the overall shape, the specific 
heat jump AC, and the ^ei{T = Q K) value. However, 
we find that the two models show qualitative differences 
for 7ei (T) at low temperatures and this low temperature 
behavior of 7ei {T) is what reflects the intrinsic properties 
of the gap symmetry, with which we can unambiguously 
identify the most compatible gap symmetry. 

A pure d-wave superconductor is expected to produce 
a T-linear specific heat coefficient such as 7(T) = aT 
because of the V-shaped density of state (DOS) N{uj) ^ 
a w. However, the large value of 7(T = 0) (~ 0.457^) in- 
dicates a substantial amount of unitary impurities. When 
the unitary impurities create a zero energy impurity band 
inside the d-wave gap, the DOS around zero energy be- 
comes a constant, N{uj) ^ Nq, which in turn produces 
a constant, temperature independent, specific heat co- 
efficient 7(r) for a finite range of low temperatures, as 
can be seen in the inset of Fig. 2 (blue circles). Our ex- 
perimental data, however, show j(T) ~ 70 -|- aT, im- 
plying that the low energy DOS should have the form 
N{uj) ^ Nq + a uj. This type of constant -|- V-shape 
DOS was indeed predicted for a ±s-wave gap with uni- 
tary impurities;^ As can be seen in Fig. 2 and its inset, 
the ±s-wave gap model with unitary impurities provide 
an excellent fit at low temperatures as well as for the 
overall shape of 7e;(r), but the d-wave model fails to 
explain the low temperature behavior. 

The dependence on magetic field of the low- 
temperature specific heat of LiFeAs is shown in Fig. 3(a) 
and 3(b) for different field orientations of H\\c and 7J_Lc, 
respectively. The lack of a Schottky-like anomaly up 
to the highest applied magnetic field enables us to un- 
ambiguously determine the Sommerfeld coefficient 'y{H) 
from the least-squares fits of C{H)/T = -i{H) + AT^ 
to the low— T specific heat. As shown in Fig. 3(c), 
the anisotropy in the quasiparticle density of states 
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FIG. 3: (color online) Low temperature specific heat capacity 
for magnetic fields applied parallel to the crystalline c-axis 
and a6-plane in panels (a) and (b), respectively. In Fig. 3(a) 
the magnetic field intensity is 0, 1, 2, 3, 5, 7, and 9 T from 
bottom to top curves, and it is 0, 1, 3, 5, 7, and 9 T in 
Fig. 3(b). Straight lines are the results from the least-squares 
fit of C/T — ^{H) + AT^ and used to obtain the zero tem- 
perature value of C/T. (c) Anisotropy of the Sommerfeld 
coefficient which is obtained from the extrapolation of the 
data in Fig. 3(a) and (b) for the two field orientations. Here 
A70 = 7o(^) ~ 7o(Or) for magnetic field along the c-axis. 



(7g(i?) - 1o)/{,1q^(H) - 70) is almost constant («2.4) 
over the experimental field range and an anisotropy ratio 
of the upper critical fields between the two field directions 
is estimated to be 1.7±0.2 through the relationship )iii^ 
(57g/(57^''=(0.5/0.3)(i?^"2Vi/^=2)°-^- The ifc2 anisotropy 
from this thermodynamic measurement is in good agree- 
ment with that obtained from transport properties where 
Hc2 was directly measured under high magnetic fieldi^i^ 

A sub- linear magnetic field dependence of 7o(i/) is dis- 
played in Fig. 4, which is often taken as evidence for a 
nodal gap structure such as a d-wave state that is known 
to have a generic \H dependence due to Doppler ef- 
fects of the nodal quasiparticles by the supercurrent cir- 
culating around the vortices. Recent theoretical work, 
however, showed that the ±s-wave state with different 
gap sizes (^smaii 7^ A/arge) cau also show a strong field 
dependence of ^{H) oc ^/~H — H?^ In Fig. 4, we show 
theoretical calculations for both ±s-wave (dotted line) 
and d-wave (dashed line) models to fit the 70 (-ff) data 
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FIG. 4: (color online) Field dependent Sommerfeld coefBcient 
extrapolated from Fig. 3a for field applied along the crys- 
talline c-axis. The black dotted line is the theoretical fit with 
the ±s-wave model and grey dashed line is the fit with the 
d-wave model, respectively. Here we used the upper critical 
field H^2 of 17 Tesla from Ref. [29]. 



with i?||c, where we use the same fitting parameters ob- 
tained from the analysis of the specific heat as a function 
of temperature (see Fig. 2). It is clear that the ±s-wave 
model equally well reproduces the field dependence of 
7(iJ) which has been taken as a signature for a nodal 
gap structure. Thus, in combination with our analysis 



of the zero field temperature dependence, our data fa- 
vors the sign changing multiband superconductivity of 
±s-wave state. This multiple SC gap scenario is con- 
sistent with recent ARPES,^^ penetration depth, ^'"^ and 
heat transport measurementSj>2£ where two distinct full 
SC gaps were reported. 

To summarize, we have reported the specific heat ca- 
pacity of a LiFeAs single crystal under magnetic field. 
The small specific heat jump ratio at Tc and substantial 
residual density of states at T = K in the SC state indi- 
cate that impurity effects due to strong unitary scatterers 
should be included to explain the specific heat results. 
The zero-field specific heat data at high temperatures 
above 2 K and sub-linear magnetic field dependence of 
the Sommerfeld coefficient are equally well described by 
nodal and nodeless multiband SC gaps with impurity ef- 
fects. However, the low-temperature electronic specific 
heat below 2 K excludes the nodal d-wave state, but is 
consistent with the sign changing multiband supercon- 
ductivity of lbs-wave state. 
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